Magnetic polarons and magnetoresistance in EuB 6 
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EuB6 is a low carrier density ferromagnet which exhibits large magnetoresistance, positive or 
negative depending on temperature. The formation of magnetic polarons just above the magnetic 
critical temperature has been suggested by spin-flip Raman scattering experiments. We find that 
the fact that EuB6 is a semimetal has to be taken into account to explain its electronic properties, 
including magnetic polarons and magnetoresistance. 



I. INTRODUCTION 

Europium hexaboride is part of the large and heteroge- 
neous class of materials that exhibit Colossal Magnetore- 
sistance (CMR). The ferromagnetic transition in EuB>6 is 
accompanied by a dramatic change in resistivity. There 
is a large body of experimental data available on the mag- 
netic and electric properties, but a thorough understand- 
ing is lacking. 

EuB>6 has a cubic unit cell with Eu-ions at its vertices 
and a Boron octahedron at its center. The material is 
ferromagnetic and shows two magnetic transitions: at 
T M = 15.3K and at T c = 12.5K Q. These have been 
associated with a spin reorientation and a ferromagnetic 
transition, respectively. Neutron diffraction experiments 
have given a magnetic moment /i ~ 7.3 + 0.5/1^. This 
is exclusively due to the localized half-filled /-shell in the 
Eu+ 2 ions. [3 

Electronic structure calculations j|, Shubnikov-de 
Haas and de Haas-van Alphen measurements Q show 
that EuBg is a semimetal. The Fermi surface consist of 
two ellipsoidal pockets, one electron-like and one hole- 
like, centered on the X point in the Brillouin zone. The 
pockets contain very few carriers: Hall effect measure- 
ments yield n c g ~ 10~ 2 carriers per formula unit Q at 
low temperatures. Small dilations of the boron octahedra 
cause overlap of the conduction and valence bands at the 
X points rendering EuBg semimetallic. The carrier con- 
centration decreases smoothly as temperature increases. 

The electrical resistivity is metallic in the ferromag- 
netic regime. It shows a sharp peak near Tq . Above this 
temperature, the resistivity decreases with an almost ac- 
tivated temperature dependence until it reaches a mini- 
mum at about 30K. At higher temperatures it increases 
and eventually starts to saturate at about room temper- 
ature. 

The application of a magnetic field produces sharp 
changes in the resistivity. Close to the magnetic tran- 
sition, negative magnetoresistance (MR) values of up to 
100% have been observed [a,l9- This decrease in resistiv- 
ity is accompanied by a large decrease in the (negative) 
Hall coefficient 5] and an increase in the plasma frequency 
• The change in the plasma frequency is more gradual 
than the changes in resistivity and Hall coefficient. 

In the ferromagnetic regime, on the other hand, the 



MR is large and positive: at 1.7 K resistivity changes of 
up to 700% have been observed in a transversal applied 
field of 7T HQ. The MR depends quadratically on the 
applied field strength at low temperatures |8|. 

Just above Tc and up to ~ 30K, the existence of mag- 
netic polarons has been proposed as the cause of the 
Stokes' shift measured with Spin Flip Raman Scatter- 
ing 

(SFRS)H[j3. 

The resistivity is activated at these 
temperatures. However, this data contains a previously 
unremarked puzzle, in that the energy scale turns out 
to be considerably lower - by a factor of thirty - than 
expected based on reliable estimates of the exchange in- 
teraction. 

Wigger et al.Jj showed how the crossover between 
large positive and large negative MR from well below to 
well above the ferromagnetic transition can be explained 
by the dominance of orbital scattering at T <C T c to spin 
scattering at T 3> Tc. The model we shall use for the car- 
rier transport in these regimes is similar to that of ref. 
and we shall thus suppress most of the details. The key 
feature of the model is its multiband nature - there are 
two types of carrier. 

In this paper we concentrate principally on the regime 
close to T c and analyze the evidence for the existence 
of magnetic polarons in Europium Hexaboride. We show 
how the SFRS results can be explained using a multiband 
model, resolving the conundrum of the anomalously small 
energy associated with the carrier spin flip. 

II. MODEL AND PARAMETERS 

We model EuBg as a ferromagnetic semimetal with a 
low carrier density. Both electrons and holes are itinerant 
and are coupled to the local moments S — 7/2. This can 
be described by the following general Hamiltonian: 

H = -t ]T c+ c> - jJ2 S t Sj 

i i 

Here, the hopping parameter is roughly t — 0.1 eV 3]. 
c ta-(e h) 1 ~<r,v' c i<r(e,h) is the itinerant carrier spin operator 
and the subindices e and h stand for electrons and holes 
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respectively. J' e (J' h ) is the on-site coupling between the 
spins of the electrons (holes) and the local moments. J 
is the magnetic exchange between local moments. 

First of all we need to discuss what is the origin of 
the ferromagnetism and the order of magnitude of the 
magnetic couplings. Ferro- and antiferro-magnetism of 
the insulating Eu-chalcogenides (EuX, X= O, S, Se, Te) 
has been explained as due to superexchange interaction 
between neighbor Eu ions [ill Il2fl through the anion 
between them. The density of carriers in the undoped 
chalcogenides is too low to expect any indirect RKKY 
(Ruderman-Kittel-Kasuya-Yosida) interaction. The fer- 
romagnetic interaction arises instead from the overlap be- 
tween the 4/- and 5<i-orbitals at different cations. This 
overlap leads to an effective exchange interaction in third 
order in perturbation theory J12j. This does not apply di- 
rectly to EuBg due to the different crystalline structure, 
but nevertheless one expects that the superexchange cou- 
pling J is small. Moreover, the increase of magnetic crit- 
ical temperature and concomitant decrease of resistivity 
under high pressures has revealed that the magnetic 
exchange in EuBg is mainly due to the RKKY interac- 
tion. Therefore, J in Eq. ^is negligible. 

The RKKY magnetic exchange is mediated by the itin- 
erant carriers via their coupling with the lattice mag- 
netic moments. An effective Heisenberg-like magnetic 
exchange can be written in terms of the local Hund's like 
exchange coupling J' [J' e or J' h ) Q 

Jeff^-^^n 2 Y,F{2k F r l ) (2) 

i 

where F(x) = - xcos ^+ slDX ; n [ s the density of carriers, 
and Ep is the Fermi energy. J e // is an oscillating func- 
tion of x but is ferromagnetic for small x. This is the 
relevant limit for Europium Hexaboride, as its low car- 
rier density implies kpr — > 0. 

To estimate the value of J' we use the mean field 
relation between Tc and J e //, T c ~ zS 2 J e ff, where 
z is the coordination number for Eu, and S is the z- 
component of the local moments. Using a critical tem- 
perature T c ~ 12 K and a parabolic approximation to 
the bands, J' ~ 0.1 eV is found, consistent with reported 
data for isolated Eu-ions|l5j. In this estimation we are 
considering that only one kind of carrier is responsible 
for the magnetism. 

III. MAGNETIC POLARONS 

When the local exchange coupling J' is large enough, 
carriers can be localized by ferromagnetic clusters and 
form composite objects called magnetic polarons. Fer- 
romagnetic polarons can exist in the low temperature 
phases of antiferromagnets but here we are interested in 
those formed in the paramagnetic phase. A necessary 
condition for the existence of magnetic polarons is that 
the density of carriers is very low compared to the in- 
verse of the correlation volume, namely n£ 3 << 1. When 



this condition is fulfilled, polarons are well-defined non- 
overlapping entities. 

There are two kinds of magnetic polarons: free and 
bound. A free magnetic polaron is a carrier localized 
in a ferromagnetic cluster embedded in a paramagnetic 
background. A carrier that is coupled strongly to local 
moments via a Hund's like coupling tends to align the 
moments that are within a Bohr radius. This causes a 
trapping potential that localizes the carrier. The poten- 
tial can be enhanced by random fluctuations of the mag- 
netization that produce an alignment of local moments 
in the carrier's vicinitv[T3.ll7j. 

The carrier thus traps itself by the magnetization it 
causes. It could increase the alignment of the local mo- 
ments and hence decrease its energy by localizing itself 
in a smaller volume. However, this would lead to an in- 
crease in kinetic energy. The quantity that determines 
the stability and size of these objects is therefore J' ft 
where J' is the coupling of the carrier spin to the local 
moments and t is the hopping parameter. The ratio J' jt 
needs to be typically larger than one to guarantee 

stability of the free magnetic polaron. 

On the other hand, in bound magnetic polarons the 
main driving force trapping the carrier is not the local 
magnetic interaction but the electrostatic potential cre- 
ated by impurities. The formation of the ferromagnetic 
cluster described above does occur. However, it is a sec- 
ond order process, as the magnitude of Hund's coupling 
is much smaller than the Coulomb interaction. 

Mean field and Monte-Carlo calculations have 
shown that magnetic polarons can exist within a temper- 
ature window above Tc whose width depends on the ra- 
tio J 1 It. At higher temperatures, magnetic fluctuations 
are strong enough to destroy the magnetic polarons. Be- 
low Tq, the condition rt£ 3 << 1 is not fulfilled and the 
polarons overlap. If a magnetic field is applied within the 
existence temperature window, the size of a polaron in- 
creases until eventually the polarons overlap and produce 
a ferromagnetic transition. 

Free and bound magnetic polarons can be differenti- 
ated by their dynamics and the resistivity they cause. 
Bound magnetic polarons are bound to an impurity in the 
system so the only way of transport is via an activated 
process: when the trapped carrier is "ionized" it is free 
to move until it is trapped by another impurity. There- 
fore they produce a resistivity p such that dp/dT < 0. 
In contrast, free magnetic polarons are able to move to 
adjacent areas when random fluctuations of the nearby 
spins produce an aligned region. There is not a barrier 
to overcome in this process. This transport mechanism 
has been called "fluctuation-induced hopping" [2(| and 
produces a metallic resistivity dp/dT > 0. 

Magnetic polarons have been largely studied in con- 
nection with Eu-chalcogenides (EuO, EuS, EuSe, EuTc) 
[2ll I22I and diluted magnetic semiconductors such 
as Cdi-^Mn^Te and Pbi-^Mn^Te with x, the concen- 
tration of magnetic ions, small. Experimental evidence 
included photoluminescence spectra [24| and magneto- 



3 



optical experiments as Spin Flip Raman Scattering [HI 
I25L l26l |. The Raman scattering spectrum shows an in- 
elastic peak at low frequencies (Stokes' shift) which, for 
the diluted magnetic semiconductors, depends as follows 
on polaronic properties: [lfj 



AE = xJ'(S z ), 



(3) 



where x is the density of magnetic ions participating in 
the formation of polarons |27|. and J' is the local ex- 
change interaction between the s itinerant electrons and 
the d electrons localized in the Mn ions. The low density 
of local moments makes for a small Stokes' shift, which - 
for Cdi_a;Mn x Se - is consistent with experiment [l6|. 

SFRS measurements done in EuB 6 have similarly re- 
vealed a zero- field peak in scattered intensity of the order 
of - \2meV UGH at 18K , J ust above the 

magnetic criti- 
cal temperature. The behavior of this peak with temper- 
ature and external magnetic field is consistent with the 
stability conditions theoretically established for magnetic 
polarons. Free magnetic polarons are not expected to be 
stable in EuBg as J' and t are comparable. Moreover, we 
have argued above that the activated behavior of the re- 
sistivity is better explained by means of bound magnetic 
polarons. Eu-site vacancies would produce the binding 
Coulomb potential for electrons. 

We expect the energy of the Stokes' shift to be given 
by Eq. [3] but now x ~ 1 as EuBg has a local moment 
on every Eu site in the cubic lattice. Very few site va- 
cancies are expected in this fairly clean material. Mean- 
field theory 22} predicts that bound magnetic polarons 
are fully spin polarised so (S z ) = 7/2. Using this value 
and the energy of the Stokes' shift we obtain J' ~ 3meV. 
This is far too low compared to the values reported in the 
literature for J' in isolated Eu ions ~ lOOmeV [9I Il5|. 

We are therefore left with a conundrum: the peak in 
the light scattering intensity follows all the trends calcu- 
lated for an object with magnetic origin but the energy 
of that peak is almost two orders of magnitude smaller 
than expected. The solution to this problem lies in the 
fact that both polarised electrons and holes are found at 
the Fermi Energy. 

Electrons and holes come from different B and Eu or- 
bitals and therefore their magnetic couplings to the lo- 
calized spin in the Eu 4/ orbitals, J' e and J' h respec- 
tively, can be very different. Electronic structure cal- 
culations Q reveal that the hole pocket comes from the 
highest intraoctahedron B 2p band. On the other hand, 
the electron pocket comes from bonding combinations of 
the cation d orbitals pointing along the cartesian direc- 
tions with some hybridization with the B atoms and some 
free-electron-like character on the (110) axes between the 
cations. In other words, the electron charge density dis- 
tribution is mainly found around the Eu ions while the 
holes are found around the B. Therefore, the coupling of 
the electrons is expected to be much larger than that of 
the holes. Consistently, Fig. 10 in Ref.0] shows a much 
larger majority-minority spin band splitting for electrons 
than for holes. 



In conclusion, we propose that the ferromagnetic or- 
dering is produced by the itinerant electrons coupled to 
the localized spin in Eu with J' e ~ lOOmeV, while the 
itinerant holes, much more weakly coupled (J' h ~ 5meV), 
account for the SFRS Stokes' shift. A corollary of this 
identification is that there is likely a much higher energy 
feature in the SFRS, so far unobserved, that correspond 
to spin-flip of the electron state. 

In the following section we will see how the existence 
of both electrons and holes is necessary to explain other 
electronic properties of EuBg. 



IV. MAGNETORESISTANCE 
A. Positive MR at T < T c 

At low temperatures the magnetoresistance is due to 
the presence of two types of carriers and we will call it 
"orbital" MR. There are two effects involved. The same 
physics that causes the Hall effect is the most important 
cause of the MR. In addition to this, there is also a small 
shift of the bands with applied magnetic field that causes 
a small change in the carrier density. 

In a simple metal with one type of carrier and a simple 
Fermi surface, there is no MR. A current that flows per- 
pendicular to a magnetic field is initially deflected due 
to the transverse Lorentz force. This causes an electric 
field, the Hall field En, to build up: 



Eu = 



-H x J. 



(4) 



The carriers that are deflected hit the edge of the sample 
and accumulate. The field that is thus built up counter- 
acts the Lorentz force. When it cancels the Lorentz force, 
the current is undeflected. This is well known; it means 
that there is no MR in a normal metal. The component 
of the current density that is parallel to the applied elec- 
tric field is not affected by the magnetic field. Therefore 
the resistivity remains unaffected as well. 

In a semimetal, on the other hand, there are by defi- 
nition two kinds of carriers. These two kinds will almost 
invariably have different scattering times and different 
masses. Eq. 0] shows that the Hall voltage depends on 
the scattering time and the mass of the carrier. There- 
fore, the Hall voltages of the different kinds of carriers 
are also different. There is thus no voltage at which the 
two carriers will travel through the sample without de- 
flection. 

The difference in the Hall voltages is proportional to 
the applied magnetic field, so that the component of the 
current parallel to the applied electric field decreases with 
increasing field. The resistivity increases therefore when 
a magnetic field is applied, and the MR is therefore posi- 
tive. We call this orbital MR, since it is due to the differ- 
ence of the masses and scattering times of the two types 
of carriers. These properties derive from particularities 
of the atomic orbitals in the material. 
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The MR can be calculated from a linearized Boltz- 
mann equation that includes a magnetic held under the 
assumption of a spherical Fermi surface. This is a stan- 
dard calculation, which can be found in [2^ and [2^. We 
quote the result for the magnetoresistance: 

A/j(H) _ cr c o- h Q e - Hh) 2 H 2 

Po (o- c + Oh) 2 + H 2 (n e a e + AihCTh) 2 ' 

where a c and <7h are the electron and hole conductivity 
respectively, and fi c and [Ah are their mobilities, both in 
the absence of any magnetic held. The right hand side of 
the above formula is easily seen to vanish if both types 
of carriers have the same mass and scattering time: the 
difference of the mobilities in the numerator vanishes in 
that case. 

At low temperatures the orbital MR is largest: the 
scattering time is largest, so that the deflection is largest 
as well. As the scattering time decreases the effect be- 
comes less important. For sufficiently small scattering 
times orbital MR becomes negligible. 

The second effect that causes MR in EuBg is the shift- 
ing of the bands when a magnetic field is applied. The 
shift is caused by the coupling J' of the carriers to the 
Eu local moments. Let us write the local moments as 
the sum of their average and the deviations therefrom: 
Si = (S)+8Si. We now use this expression in the Hamil- 
tonian in Eq. ^ and obtain a term that couples the car- 
rier's energy to the magnetization. As the magnetization 
grows, the electron-like band is shifted to lower energies, 
and the hole-like band to higher energies. This causes an 
increase in the number of carriers as carriers spill over 
from one band into the other. 

The change in carrier density can be obtained from the 
following two requirements. Firstly, the band-shift intro- 
duced by the change in magnetization M is 8E = J' e M/2 
for electrons and 8E = J' h M/2 for holes. Secondly, the 
increase in the number of carriers of both kinds is equal: 
as the hole like band is shifted up, new holes are created 
as negatively charged particles spill into the electron-like 
band and vice versa. Overall charge balance is main- 
tained. We also assume a spherical Fermi surface, an 
accurate assumption in the case of EuBg [l| . We can use 
the requirement that charge neutrality be conserved to 
calculate the shift in the Fermi level. The shift is then 
used to obtain the number of carriers by integrating the 
density of states up to the Fermi level. The increase re- 
sulting from the shift of the bands is small, even at full 
saturation of the magnetization. Its effect on the MR at 
low temperatures is then negligible. The magnetization 
as a function of applied field is obtained from a Curie- 
Weiss model. 

We included this change in carrier density due to the 
bands shifting in our model for the orbital MR. The con- 
tribution to the MR due to band shifting is opposite to 
that of the orbital effects. An external magnetic field in- 
creases the carrier density, so that it decreases the resis- 
tivity. Since the magnetization is almost saturated at low 
temperatures, the carrier density does not change much 
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FIG. 1: Orbital MR at low temperatures for different applied 
magnetic fields. The data is taken from ref[f|. The plateau 
seen in the model below 5 K is an artifact of data digitization. 



with applied field, and the MR is affected only slightly. 
The orbital contribution to the MR dominates. 

We plot the MR obtained for the combined effect of 
the band shifts and the orbital effects in Fig. ^ The 
temperature dependent scattering times at zero field are 
obtained from experimental data in [3 • The electron and 
hole mobilities are obtained from the conductivity at zero 
field, the masses of the carriers Q and the carrier densi- 
ties. The latter were obtained from the plasma frequency 
in [3 . We introduced a small imbalance between the car- 
rier densities of 6.10 _4 g e per unit cell, in accordance with 
[4j. This imbalance is thought to arise from impurities. 
These numbers provide input to the model. 

Our simple model, which depends only on parameters 
measured at zero field, can reproduce the large positive 
MR at low temperatures accurately. Fig. E] shows that it 
reproduces the magnitude of the MR well. It also predicts 
the (nearly) quadratic dependence in applied held. From 
Fig. it is clear that the orbital contribution to the MR 
dwindles at higher temperatures. The model proposed 
will lose its validity near the magnetic transition and in 
the paramagnetic phase, when other effects dominate. 



B. Negative MR at T » T c and above 

The mechanisms that govern the low-temperature MR 
become insignificant near the ferromagnetic transition. 
Close to T c the scattering time is so short that the pos- 
itive orbital contribution to the MR is negligible. On 
the other hand, the shift of the bands caused by an ap- 
plied magnetic field becomes substantial. We estimate 
that the carrier density changes by about 7% as the ap- 
plied field saturates the magnetization. This could ex- 
plain only part of the negative MR in the critical regime. 

Aditionally, near the critical point, spin fluctuations 
may provide a large contribution to the electrical resis- 
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tance. The dominant modes near a ferromagnetic transi- 
tion are those with small q. They only produce substan- 
tial backscattering if 2kF is itself small • This is 
the case of EuB6 as its carrier density is very small. The 
suppression of the spin fluctuations when a magnetic field 
is applied is largely responsible for the MR found in the 
critical regime. Moreover, as shown in Ref. 01 ( see Fig- 
3), the localization of the carriers in magnetic polarons 
further increases the magnetoresistance. 

In the temperature regime - just above T c - where the 
SFRS data gives evidence for bound magnetic polarons, 
the magnetoresistivity is large, and strongly negative, 
as expected since an applied magnetic field suppressed 
the magnetic polarons. At temperatures greater than 
about 30K where the polarons are destabilized, we have 
a smaller (but still negative) MR dominated by local spin 
fluctuation scattering. 

As mentioned in the introduction, the negative MR 
is accompanied by a shift in the plasma frequency 0. 
The carrier density change produced by the band shift- 
ing alone is too small to account for that shift. This is 
in contrast to the results by Wigger et al [8j who claim 
consistency between the experiment and their calcula- 
tion, though the functional dependence M 2 oc (tu 2 ) 2 
is not reproduced. There are two main differences be- 
tween their model and ours: they consider EuBg to be 
a strongly compensated n-type magnetic semiconductor 
and use the same local coupling J' for electrons and holes. 



V. SUMMARY 

EuBg is a low carrier density ferromagnet with unusual 
properties: the resistance changes from metallic to acti- 
vated and then to metallic as the temperature increases 
and the magnetoresistance changes sign close to Tc- The 
activated region has been ascribed to the existence of 
bound magnetic polarons. We discuss their existence in 
the light of Spin Flip Raman Scattering measurements re- 
ported in Ref. [UEI . We conclude that the signature seen 
by those experiments is due to the Hund's like coupling of 
holes with the local spins while electrons are responsible 
for the magnetic ordering through the RKKY interaction. 
This resolves the puzzle that the RKKY transition tem- 
perature implies an exchange coupling of the carriers to 
the local moment of about 0.1 eV, about 30 times larger 
than the measured spin-flip energy of a carrier trapped 
in a bound polaron. The existence of bound magnetic 
polarons is also consistent with a large negative magne- 
toresistance above T c . The positive magnetoresistance in 
the ferromagnetic phase is also produced by the inter- 
play of two kinds of carriers with different masses and 
scattering times. 
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